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Using the effective potential approach for composite operators we have formulated the quantum 
model of the QCD vacuum. It is based on the existence and importance of the nonperturbative 
g~ 4 -type dynamical, topologically nontrivial excitations of the gluon field configurations. The QCD 
vacuum is found stable since the vacuum energy density has no imaginary part. Moreover, a possible 
ground (stationary) state of the nonperturbative Yang-Mills (quenched QCD) vacuum is discovered. 
The vacuum energy density at stationary state depends on a scale at which nonperturbative effects 
become important. The quark part of the vacuum energy density depends in addition on the constant 
of integration of the corresponding Schwinger-Dyson equation. The value of the above mentioned 
scale is determined from the bounds for the pion decay constant in the chiral limit. Our value for 
the chiral QCD vacuum energy density is one order of magnitude bigger than the instanton based 
models can provide while a fair agreement with recent phenomenological and lattice results for the 
, chiral condensate is obtained. 
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I. INTRODUCTION 



The nonperturbative QCD vacuum has a very rich dynamical and topological structure. It is a very complicated 
medium and its dynamical and topological complexity means that its structure can be organized at various levels 
(classical, quantum) and it can contain many different components and ingredients which contribute to the vacuum 
energy density, the one of main characteristics of the QCD ground state. Many models of the QCD vacuum involve 
some extra classical color field configurations such as randomly oriented domains of constant color magnetic fields [1], 
background gauge fields, averaged over spin and color [2], stochastic colored background fields [3], etc (see also Ref. 
[4] and references therein). The most elaborate random and interacting instanton liquid models (RILM and IILM) 
of the QCD vacuum [5] is based on the existence of the topologically nontrivial instanton-type fluctuations of gluon 
fields, which are solutions to the classical equations of motion in Euclidean space [6-8]. 

Today there are no doubts left that the dynamical mechanisms of the important nonperturbative quantum phenom- 
ena such as quark confinement and dynamical (or equivalently spontaneous) chiral symmetry breaking (DCSB) are 
closely related to the complicated topologically nontrivial structure of the QCD vacuum [9-12]. On the other hand, 
it also becomes clear that the nonperturbative infrared (IR) dynamical singularities, closely related to the nontrivial 
vacuum structure, play an important role in the large distance behaviour of QCD [13,14]. For this reason, any correct 
nonperturbative model of quark confinement and DCSB necessarily turns out to be a model of the true QCD vacuum 
and the other way around. Our model of the true QCD ground state is based on the existence and importance of 
such kind of the nonperturbative, quantum excitations of the gluon field configurations which effectevely correctly 
can be described by the q~ 4 behaviour of the full gluon propagator in the IR domain (at small q 2 ). It describes the 
zero modes enhancement (ZME) dynamical effect in QCD at large distances [15-19] (for additional references see Ref. 
[19]). These excitations are also topologically nontrivial in comparison with the free gluon structure, q~ 2 ■ 

In this context, we note that the attractive classical model of the QCD vacuum as a condensation of the color- 
magnetic monopoles (QCD vacuum is a chromomagnetic superconductor) proposed by Nambu, Mandelstam and 't 
Hooft and developed by Nair and Roscnzweig (see Ref. [20] and references therein) as well as the classical mechanism 
of the confining medium [21] and an effective theory for the QCD vacuum proposed in [22], also invoke q~ 4 behaviour 
of the gluon fields in the IR. Let us underline that without q~ 4 component in the decomposition of the full gluon 
propagator it is impossible to obtain the area law for static quarks (indicative of confinement) within the Wilson 
loop approach [23] . This behaviour of the full gluon propagator in the IR is also required to derive the heavy quark 
potential within the recently proposed exact renormalization group approach [24]. 

In the next section for the sake of the reader's convenience we provide a brief review on the dynamical equations 
approach to nonperturbative QCD, since our quantum model of the QCD vacuum is its direct consequence. 
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II. DYNAMICAL EQUATIONS APPROACH TO NONPERTURBATIVE QCD 



Our approach to nonperturbative QCD is based on solutions to the quark and ghost Schwinger-Dyson (SD) quantum 
equations of motion which should be complemented by the investigation of the corresponding Slavnov- Taylor (ST) 
identities [19,25,26]. If it is true that QCD is an IR unstable theory (i. e. if it has no IR stable fixed point) [27] then 
the low-frequency modes of the Yang-Mills (YM) fields might be enhanced due to nonperturbative IR divergences 
[28]. So the full gluon propagator can diverge faster than the free one at small momenta [15-29] 

DM ~ (<7 2 r 2 , q 2 ^0, (2.1) 

which might be equivalently refered to as the strong coupling regime [27]. If indeed the low- frequency components of 
the virtual fields in the true vacuum have a larger amplitude than those of the bare (perturbative) vacuum [15], then 
the Green function, describining the propagation of a single quark in the true QCD vacuum, should be reconstucted 
on the basis of this effect. It is important to understand that a possible ZME effect (2.1) is our primery dynamical 
assumption.^ We consider this effect as a very similar confining ansatz for the full gluon propagator in order to use it 
as input information for the quark and ghost SD equations as well as in the corresponding quark-gluon ST identitiy. 

Such a singular behaviour of the full gluon propagator in the IR region requires the introduction of a small IR 
regularization parameter e, in order to define the initial equations in the IR by the dimensional regularization method 
[30] within the distribution theory [31]. This yeilds the regularization expansion for the above mentioned strong IR 
singularity as follows (of course, under integrals in four dimensional Euclidean space) [19,25,26,31] 

(q 2 )~ 2+c = —5 4 {q) + finite terms, e -> 0+ (2.2) 

and the terms of order e are not shown here for simplicity. Because of this, the quark propagator and other Green's 
functions become dependent, in general, on this IR regulation parameter e, which is to be set to zero at the end of 
computations, e — > + . For the sake of brevity, this dependence is always understood but not indicated explicitly. 

There are only two different types of behaviour of the quark propagator with respect to e in the e — > + limit. If 
the quark propagator does not depend on the e - parameter in the e — ► + limit then one obtains the IR regularized 
(from the very beginning) quark propagator. In this case quark confinement is understood as the disappearance of 
the quark propagator pole on the real axis at the point p 2 — m 2 , where m is the quark mass. Such an understanding 
(interpretation) of quark confinement comes, apparently from Preparata's massive quark model (MQM) [32] in which 
external quark legs were approximated by entire functions. A quark propagator may or may not be an entire function, 
but in any case the pole of the first order (like the electron Green's function has in QED) disappears (see Ref. [19] 
and references therein). However, the absence of the pole — type singularities in the quark Green's functions is 
only the first necessary condition of the quark confinement at the fundamental (microscopic) quark-gluon level. At 
hadron (macroscopic) level there exists the second sufficient condition, namely the corresponding Bethe-Salpeter 
(BS) equation for the bound-states should have the discrete spectrum only [33] in order to prevent quarks to appear 
in asymptotically free states. Consequently, the meson-quark decay amplitude vanishes in the IR limit e — > + , 
consistent with confinement [18]. Thus in general case the confinement criterion consists of the two above formulated 
conditions. This definition generalises the linearly rising potential between heavy quarks since it is relevant not only 
for light quarks but for heavy quarks as well. On the other hand, a quark propagator can vanish after the removal 
(e — > + ) of the IR regularization parameter e. A vanishing quark propagator is also a direct manifestation of quark 
confinement (see again Ref. [19] and references therein). Such understanding of quark confinement comes, apparently, 
from two-dimensional QCD with N c large limit [33]. 

Here a few remarks are made in advance. From our approach to nonperturbative QCD, it follows that the first 
nontrivial approximation to the quark-gluon vertex (as it is required by the correct treatment (2.2) of the above 
mentioned strong IR singularity (2.1) by the distribution theory [31] within the dimensional regularization method 
[30]) is the quark-gluon vertex at zero momentum transfer. It depends on four independent form factors which should 
be determined within the corresponding ST identity. Obviously, in this case there is no sense to speak about potential 
at all (in particular linearly rising one between heavy quarks) and it is necessary only to speak about the large scale 



1 However, let us remind the reader that after the pioneering papers of Mandelstam in the covariant (Landau) gauge [15] and 
Baker, Ball and Zachariasen in the axial gauge [16], the consistancy of the singular asymptotics (2.1) with direct solution to 
the SD equation for the full gluon propagator in the IR domain was reapeatedly confirmed (see for example Refs. [17,19,29] 
and references therein). 
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structure of the QCD vacuum and its revelance to quark confinement and other nonperturbative effects in QCD. fn 
other words, nobody can believe that the potential concept is useful in this case.0 This vertex might be approximated 
by point-like vertex only for heavy quarks when one neglects corrections to the vertex induced by virtual gluons. In 
this case the linearly rising potential becomes relevant and that is why it was "seen" in lattice QCD [34]. This should 
be considered as strong lattice evidence (though not direct) of the existence and importance of g _4 -type excitations 
of gluon field configurations in the QCD vacuum. As it has been already mentioned above, without this component in 
the decomposition of the full gluon propagator in continuum theory it is impossible to "see" linearly rising potential 
by lattice simulations not involving some extra (besides gluons and quarks) degrees of freedom. 

We develop a method for the extraction of the IR finite (regularized) Green's functions in QCD. The IR finiteness 
of the Green's functions means that they exist as e — > + . For this purpose, we have worked out a renormalization 
program in order to cancel all IR nonperturbative divergences, which makes it possible to explicitly show that all 
Green's functions are IR multiplicative renormalizable (IRMR). We have also shown that for the covariant gauges 
the complications due to ghost contributions can be considered in our approach. A closed set of equations for the IR 
finite from the very beginning quark propagator (S = S) and other IR renormalizcd quantities in the quark sector is 
worked out to be [19,25,26] 

S-\ P ) = S^(p)+fi 2 r^p,0)S( P h, l + S uv (p), (2.3) 
r p (p, 0) = idpS- 1 (p) - S(p)T„ (p, 0)5- 1 (p) + T u / (p, 0) (2.4) 

with d^ = d/dp^. Here So(p) and S(p) are the free and full quark propagators, respectively while r M (p, 0) is the 
corresponding quark-gluon full vertex at zero momentum transfer. Evidently, the second terms in Eqs. (2.3-2.4) are 
due to the deep IR (q~ 4 ), while the third terms are due to the ultraviolet (UV) components in the decomposition of 
the full gluon propagator. The characteristic IR renormalized mass scale parameter, fi 2 , determines in general the 
physical scale of nonperturbative dynamics within our approach. In particular, it is directly responsible for dynamical 
breakdown of chiral symmetry at the fundamental quark level.^j Indeed, if it is zero (/2 2 = 0) then from Eq. (2.3) 
it follows that there is no contribution to the quark mass generated function (quark self-energy) from the deep IR 
(confinement) region. In other words, the nonperturbative phase does not exist at all in this case. It comes from the 
initial mass scale parameter u 2 which is necessary to introduce due to the q~ A behaviour of the full gluon propagator 
in the deep IR region. The IR finite quark wave function renormalization constant which should multiply the free 
quark propagator in Eq.(2.3) is to be set to one without loosing generality because of the above mentioned IRMR 
property of our approach [19,25]. 

We first approximate the exact quark SD equation (2.3-2.4) by its deep IR (confinement) piece only assuming that 
precisely this term is mainly responsible for nonperturbative effects in QCD in particular quark confinement. That is 
why we do not present the explicit, rather complicated expressions for the corresponding UV terms in Eq. (2.3-2.4). 
Because of the same reason, the quark wave function UVMR constant is also set to one since the confinement piece 
of the quark SD equation is free from UV divergences. The nontrivial UVMR program should be done either from 
the very beginning or after completing our IRMR program when one takes into account the UV terms as well. Let us 
note in advance that approximating the full gluon propagator by its deep IR (confinement) piece only in the whole 
range, we nevertheless will obtain a solution for the dynamically generated quark mass function, which manifests the 
existence of the effective boundary value momentum (sec below). 

Introducing then the appropriate Euclidean dimensionless momentum variable x — p 2 / p 2 , and doing some algebra, 
the quark SD equation (2.3-2.4) (without UV pieces) in terms of the quark propagator iS(p) = pA(p 2 ) — B(p 2 ) form 
factors, becomes [19,25,26] 

xA' = -(2 + x)A-l-m B, (2.5) 
2BB' = -3A 2 - 2{B - m A)B, (2.6) 

where A, B = A{x), B{x) and the prime denotes the derivative with respect to x and p = 7^^. Let us emphasize 
that the obtained system of equations (2.5-2.6) certainly leads to dynamical breakdown of chiral symmetry (quark 
mass generation) since a chiral symmetry violating solution (mo = 0, A{x) ^ 0, B(x) ^ 0) is only allowed, while 
a chiral symmetry preserving solution (mo = B(x) = 0, A(x) ^ 0) is forbidden. 



2 One of the authors (V.G.) is grateful to G.'t Hooft for the discussion of the properties of the q~ 4 singularities from this point 
of view. 

3 Precisely this interpretation will be put into the heart of our scale-setting scheme in order to perform numerical calculations, 
see section 4 below. 
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It is easy to check that solutions to this system in the chiral limit (m = 0) are 

A(x) =x- 2 {l-x-e- x ), (2.7) 
B 2 (x , x) = ie- 2x f X ° dx'e 2x 'A 2 (x'), (2.8) 

J X 

where xo = Pq/JI 2 is the corresponding (arbitrary at this stage) constant of integration. It plays the role of UV cut-off 
(dimensionless) while po is the UV cut-off in momentum space. It is easy to see that obtained solution for the quark 
propagator is regular at zero, has no pole- type singularities (indicative of confinement). This remains valid for the 
nonchiral case as well. As it has been already underlined above, it corresponds to dynamical breakdown of chiral 
symmetry (quark mass generation). It is also nonperturbative ( it can not be expanded in powers of the coupling 
constant) and the function A(x) automatically approaches the free propagator at infinity (asymptotic freedom). In 
the solution of the differential equation for it, the constant of integration was put to infinity from the very beginning 
unlike for the quark mass generated function B(x). 

The remarkable feature of the solution (2.8) for the dynamically generated quark mass function is that it exhibits 
an algebraic branch points at x — xq and at infinity (x — > oo at fixed xq), which are caused by the inevitable ghost 
contributions in the covarint gauge. That is why the limit x — > is unphysical since the dynamically generated quark 
mass function J3 2 (x ,x) becomes pure imaginary which contradicts the idea that the confining particle should have 
no imaginary part [35]. On the other hand, this limit is achieved when p, 2 — > oo (at fixed po), but this is impossible 
since the physical scale of the nonperturbative dynamics jl is either finite (then the nonperturbative phase exists) 
or zero (only the perturbative phase survives, see Eq. (2.3)), so it can not be arbitrarily large. The perturbative 
(UV) limit xo — > oo in terms of the constant of integration xo = Po/fi 2 can be achieved by two different ways. First, 
it is recovered at fixed p 2 , when the renormalized mass scale parameter p, 2 goes to zero (/I 2 — ► 0) indicating that in 
this limit there is no nonperturbative phase at all and only the perturbative phase remains, see Eq. (2.3). Second, 
the UV limit is also recovered at fixed jl 2 , but the UV cut-off p n goes to infinity (p — > oo). However, one can not 
put automatically x = oo for the dynamically generated quark mass function B 2 (xo, x) since the integral (2.8) with 
xq = oo for any finite x in particular x — does not exist at all. Hence, we have to keep the constant of integration Xo 
in (2.8) arbitrary but finite in order to obtain a regular and nontrivial solution for the deep IR region, x > x. Thus 
it plays the role of the UV cut-off, separating deep IR (confinement) region, where nonperturbative effects become 
dominant, from the intermadiate and UV regions where they might be neglected. From the solution (2.8) it also 
follows that extrapolating it to infinity (x — > oo), it should be simultanously accompanied by the UV limit in terms of 
the constant of integration xo — > oo in order to avoid the influence of the unphysical singularity (the above mentioned 
branch point at infinity at fixed xo) on our solutions. In other words, the correct UV limit x — > oo requires x — > oo 
and vice versa. In this case the dynamically generated quark mass function (2.8) identically vanishes in accordance 
with the vanishing current quark mass in the chiral limit. In what follows in all integrals containing the quark degrees 
of freedom (the dynamically generated quark mass function i3 2 (x ,x)) explicitly, the perturbative (UV) limit shoud 
be understood in this sense. Thus one may conclude, that the analytical properties of our solutions (structure of the 
singularities) influence both the IR and UV behaviour of the quark propagator. 

Concluding this section, we note that our solution does not explicitly depend on ghost degrees of freedom as well 
as on gauge choice. It is worthwile to emphasize in advance, that within our approach [19,25,26] the explicit ghost 
and gauge dependence is shifted from confinement piece to the UV piece of the full quark propagator (2.3-2.4). 

III. ZERO MODES ENHANCEMENT QUANTUM MODEL OF THE QCD VACUUM 

As it was mentioned in Introduction any correct nonperturbative model of quark confinement and DCSB necessarily 
becomes a model of the QCD ground state, i. e. its nonperturbative vacuum. The effective potential approach for 
composite operators [36,37] allows us to investigate the QCD vacuum, since in the absence of external sources the 
effective potential is nothing but the vacuum energy density, one of the main characteristics of the nonperturbative 
vacuum. It gives the vacuum energy density in the form of loops expansion where the number of 2PI vacuum loop (con- 
sisting of confining quarks with dynamically generated quark masses and nonperturbative gluons properly regularized 
with the help of ghost) is equal to the power of the Plank constant, h. So here we will establish completely quantum 
part of the vacuum energy density which is due to g~ 4 -type nonperturbative, topologically nontrivial excitations of 
the gluon field configurations in the QCD vacuum. 

It is convenient to start from the quark part of the vacuum energy density which to leading order (log-loop level 
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~ fi-)0is given by the effective potential for composite operators as follows [36] 

j^-Tr {HS^S) (S^S) + 1} , (3.1) 

where S(p) and So(p) are the full and free quark propagators, respectively. Here and everywhere below in this section 
the traces over space-time and color group indices are understood. Let us note that the effective potential (3.1) is 
normalized as V(Sq) = 0, i. e. the perturbative vacuum is normalized to zero. In order to evaluate the effective 
potential (3.1) we use the well-known expression, 

TrlniS^S) = 3 x bxdet^S) = 3 x 2 hip 2 [p 2 A 2 {-p 2 ) - B 2 {-p 2 )] , (3.2) 

where p 2 A 2 (-p 2 ) - B 2 (-p 2 ) = y/det[-iS(p)]. The factor 3 conies from the trace over quark color indices. Going over 
to Euclidean space (n — 4, d A p — > i^p, p 2 — » — p 2 ), in terms of dimensionless variables and functions (2.5-2.8), we 
finally obtain after some algebra (e q = V(A, £>)), 

1 3 

= — e i = FT x x a 2l g( x o,0), (3.3) 



and 



l q (x ,0)= I dxx {ln(x[xA 2 (x) + B 2 (x ,x)}) + 2xA(x) + 2} , (3.4) 
o 

where for further aims it is convenient to introduce the quark effective potential at fixed po, which becomes a function 
of dimensionless variable Xq only. Here we need to identify the UV cut-off with the constant of integration xq in 
order to guarantee that the unphysical singularities (algebraic branch points, mentioned above in section 2) in the 
B 2 (xq,x) function (2.8) will not affect the effective potential, which should be always real in order to avoid the 
vacuum instability [38] . The constant of integration xq is related to the scale below which nonperturbative effects 
become essential. For this reason, within our approach to QCD at large distances in order to obtain numerical values 
of any physical quantity, e.g. the pion decay constant (see below and Ref. [39]), the integration over the whole range 
[0, oo) reduces to the integration over the nonperturbative region [0, xo], which determines the range of validity of 
the corresponding solutions (2.7) and (2.8) for the confinement piece of the full quark propagator. We emphasize that 
the main contribution to the values of the physical quantities comes from the nonperturbative region (large distances), 
whereas the contributions from the short and intermediate distances (perturbative region), because of less singular 
behaviour in the IR, can only be treated as perturbative corrections. We shall confirm this physically reasonable 
assertion numerically. 

In order to correctly recover the perturbative (UV) limit, it is necessary to neglect the dynamically generated 
quark mass function B 2 (xo,x) (2.8) in the integrand function (3.4). As it was explained above (section 2) the correct 
perturbative limit xq — > oo requires the simultaneous UV limit, x — + oo. Then one obtains that fi g as a function 
of xq at upper limit (xq — » oo) converges as Xq lnxo, satisfying thereby the initial normalization condition of the 
perturbative vacuum to be zero. At the same time, the quark effective potential (3.3-3.4) as a function of xq has no 
stationary state (a local minimum) while the regularized gluon part does have (see below). Let us note also that at 
fixed p,, the corresponding effective potential Cl q = (1//2 )e q as a function of xq has uncorrect UV limit (it diverges as 
lnrro at upper limit, xq — > oo) contradicting thereby the above mentioned normalization condition. Thus ft can not 
be fixed. It should be determined from some good physical observable (see section 4 below). The difference in the 
behaviour of the effective potential (as a function of xq) at fixed po or p, should be traced back to its two versions 
when the variable is the quark propagator, S, [36] or the quark self-energy, S, [37], respectively. 

We evaluate the nonperturbative gluon part of the effective potential, which at the log-loop level is given by [36] 



V{D) = \\ Wr Tr{HD ° lD) {D ° lD) + 1} ' 



(3.5) 



4 Next-to-leading and higher terms (two and more vacuum loops) are suppressed by one order of magnitude in powers of h at 
least and are left for consideration elsewhere. 



5 



where D(p) is the full gluon propagator and Dq(p) is its free (perturbative) counterpart. The effective potential is 
normalized as V(Dq) = 0, i. e. as in the quark case the perturbative vacuum is normalized to zero. In a similar way 
to Eq. (3.2), we obtain 



Tr ln(L> 1 L>) = 8 x lndet(D 1 D) = 8 x 4 In 



I^ 2 ) + \ 



(3.6) 



where the factor 8 is due to the trace over the gluon colour indices and it becomes zero (in accordance with the above 
mentioned normalization condition) when the full gluon form factor is replaced by its free counterpart by setting 
simply d(—p 2 ,a) = 1 in (3.6). Let us note that this decomposition does not explicitly depend on a gauge choice. 
Approximating now the full gluon form factor by its deep IR (confinement) piece, namely d(—p 2 ) = p 2 /(—p 2 ) and 
after doing some algebra in terms of new variables and parameters (2.5-2.8), we finally obtain (in Euclidean space, 
n = 4) the following expression for the vacuum energy density due to nonperturbative gluon contributions, e g = V(D), 

tg = ^oV* x^o,0), (3.7) 

7T 



where 



zo 

I«(z , 0) = J dxx{lnx ln(x + 3) + J- - a} = | ln(l + |) - \z, \z 2 ln(l + 1) - % 2 (3.8) 
o 

and a — (3/4) — 21n2 = —0.6363. In this equation zq is, of course, the corresponding UV cut-off which in general 
differs (i. e. independent) from the quark cut-off xq, i. e. zo = q 2 /p, 2 . 

The effective potential at the log-loop level for the ghost degrees of freedom is [36] 

^2_Tr{\n{G^G) (G^G) + 1}, (3.9) 

where G(p) is the full ghost propagator and Go(p) is its free (perturbative) counterparts. The effective potential V(G) 
is normalized as V{Gq) — 0, i.e. here like in quark and gluon cases the energy of the perturbative vacuum is set to 
zero. Evaluating the ghost term e g h = V[G) in (3.9) in a very similar way, we obtain e g h = r K~ 2 k\y a ~ 2 x I g h(yo,0), 
where the integral I g h (yo , 0) depends on the IR renormalized ghost self-energy, which remains arbitrary (unknown) 
within our approach. We have introduced the ghost UV cut-off y as well, yo = k 2 /fl 2 . 

In principle, we must sum up all contributions in order to obtain total vacuum energy density. However, e g and 
e g h are divergent and therefore they depend completely on arbitrary UV cut-offs zo and yo which are in general 
different from x . We know how to calculate x by relating it to some good physical observable (see section 4 below). 
At this stage we have no idea how to calculate Zq (see however below). That is why the only way to calculate zq 
is to identify it with xq by using the arbitrariness of ghost degrees of freedom. Thus the sum e g + e g h should be 
regularized in order to define finite e g e9 , which will depend only on xo and at the same time it should be a vanishing 
function of xo at xo — > oo because of the normalization of the perturbative vacuum to zero. Approximating the full 
gluon propagator by its deep IR (confinement) asymptotics only in (3.6), the nonzero constant a appears in (3.8) 
which precisely violates the above mentioned normalization condition of the perturbative vacuum to be zero. So it 
should be additionally subtracted. For this purpose, we decompose the integral (3.8) into the three parts as follows 
l g (z ,0) — Ig (x ,0) + I g (zQ,x ) = I g (x o ,0) - (a/2)xl + Ig(z ,x Q ). Using the arbitrariness of the ghost term (3.9), let 
us subtract the unknown integral — (a/2)xQ + Ig(zo, xq) from the gluon part of the vacuum energy density by imposing 

the following condition A = — (o/2)xq+ I g (zo, x o) + Igh(yo, 0) = (having in mind that p, 4 = PqXq 2 = QqZq 2 = k^y^ 2 ). 
In fact, we regularize the gluon contribution to the vacuum energy density by subtracting unknown term by means 
of another unknown (arbitrary) ghost contribution, i. e. e g + e g h — e r g eg — tt~ 2 PqXq 2 I g (xo, 0) and in what follows 
superscript "reg" will be omitted for simplicity. It is easy to show that the above mentioned condition of cancellation 
A = of the UV divergences (due to arbitrary zo, yo) is consistent with this definition of the regularized vacuum 
energy density of the nonperturbative gluon part. At the same time, it becomes negative and has a local minimum 
(see below), while the unregularized contribution tt~ 2 PqXq 2 I g (xo, 0) is always positive (because of constant a). Thus 
our rcgularization procedure is in agreement with general physical interpretation of ghosts to cancel the effects of the 
unphysical degrees of freedom of the gauge bosons [27,40]. 

The regularized vacuum energy density due to the nonperturbative gluon contributions becomes 

" 9 = pX 9 = i X x °~ 2 {91n(l + y) - l*o - 4ln(l + A)}, (3.10) 
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where, in a similar way as for quarks, we introduce the gluon effective potential at fixed poi namely tt g . Its behaviour 
as a function of xq is shown in Fig. 1. ft has a few remarkable features. First, a local minimum appears atQ 

£™ m =41n(f + -2— ) = 2.20 (3.11) 

o 

(stationary condition) and it has zero at x[j = 0.725. Second it asymptotically vanishes as xq — > oo in agreement with 
the normalization condition.^ As it was discussed above (section 2) the opposite limit xq — > (/I 2 — > oo) is unphysical 
since the nonperturbatibe IR renormalized scale p is zero or finite, it can not be arbitrarily large. That is why the 
vacuum energy density in this limit becomes positive (see Fig. 1). This obviously means that the physical region for 
parameter xq is bounded from below, namely xq > Xq = 0.725 in complete agreement with the discussion in section 
2. In this region the effective potential due to nonperturbative gluons is always negative. Thus it has no imaginary 
part (our vacuum is always stable) and at stationary state it is ri ff (x™ m ) = —0.0263 (see Fig.l). Let us also undclinc 
that without the existence of stationary state in the YM vacuum it would be impossible to calculate this number. It 
becomes clear now that this minimum would not be changed if the effective potential = (l/qfyeg would remain 
a function of zo in (3.7). What is important indeed, is to subtract the constant a from (3.8) in order to proceed to 
(3.10) as a function of Zq instead of Xo- 

The vacuum energy density for quenced QCD due to nonperturbative gluons at stationary state becomes 

t 9 = tg(xo"\Po) = -0.02634 (3.12) 

At first sight it is badly divergent since po is formally the UV cut-off in momentum space. However, within our model 
it is effectively determines the range of validity of the corresponding q~* behaviour of the full gluon propagator in 
the deep IR region (po = <7o) which, obviously can not be arbitrary large. The finiteness of this scale means that it 
might be interpreted as the effective (possibly confinement) scale of nonperturbative dynamics within our approach. 
If QCD confines such an effective scale should certainly exist. It is clear that it can not be numerically determined 
within the YM (quenced QCD) theory alone. However, it seems to us that 1 GeV is a realistic upper bound for the 
effective scale responsible for nonperturbative dynamics in quenced QCD. Let us underline, that its value has nothing 
to do with the values chosen to analyze numerical results in phenomenology or lattice approach (for example, 1 GeV 
or 2 GeV) which have no physical sense and are simply convention, while our scale has a direct and clear physical 
meaning as separating in general the nonperturbative phase from the perturbative one. 

The total vacuum energy density to leading order (as it is determined by ZME quantum model of the QCD vacuum) 
finally becomes 

e = e g + N f e q = -0.1273/2 4 + N f -^p 4 x I q (x , 0), (3.13) 

where l q (xo,0) is given in (3.4) and Nf is a number of light flavors. In derivation of the YM part of the vacuum 
energy density, e g , we use the relation p\ = x™ m p 2 (which once more shows the finiteness of po since p is always finite) 
in order to express now the vacuum energy density in terms of the fundamental scale parameter p and the constant 
of integration xq which in general need not to be numerically the same as in quenched QCD (x™ m ). That is why 
the nonperturbative effective scales for quenched (Nf — 0) and full QCD are different while the fundamental scale 
parameter of our approach p is unique and its numerical value is to be found from good physical observable in full 
QCD only. The corresponding total effective potential at fixed po, Q = — 0.1273xq 2 + NfQ q , has no local minimum. 
Thus the injection of quark degrees of freedom lifts the nonperturbative quenched QCD vacuum from its stationarity. 
How to choose the scale-setting scheme in order to numerically determine it is the subject for the next section. 

IV. THE SCALE-SETTING SCHEME 

The main problem now is to set the scale at which our calculations should be done. In our approach it means 
the choice of a reasonable value for the nonperturbative scale parameter p in (3.13). Only the requirement is that a 



5 In our previous publication, Ref. [26], the existence of the stationary state in pure gluodynamics (quenched QCD) was not, 
unfortunately, noticed. Though the quark part of the vacuum energy density was estimated correctly, the nonperturbative 
gluon contribution was, for the above mentioned reason, substantially underestimated. Stationarity of the nonperturbative YM 
vacuum was first observed and used for preliminary numerical calculations in Ref. [41]. 

6 Simiral to quark case, the gluon effective potential Cl g = (l/p 4 )e g at fixed p, has unccorect UV behaviour (it diverges as 
~ — xq as xo oo). 
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scale-setting scheme should be physically well-motiviated since ft determines the physical scale of the nonpcrturbative 
dynamics in our approach, and it should not be arbitrarily large. In our previous publications [19,26,39] the expressions 
for basic chiral QCD parameters such as pion decay constant, F (= F°), the dynamically generated quark mass 
(defined as the inverse of the quark propagator at zero momentum) and quark condensate (qq) were derived: 

3 ,2 . 3 _ 2 f Xo xB 2 (x Q ,x) 



, -^MW^M^ ^pgTBM' (41) 

m d = fi{B 2 (x ,0)y 1/2 , (4.2) 

3 3 f x ° 

{qq)o = ~^2^ x J con (a;o,0) = ~^M 3 x J dxxB(x ,x). (4.3) 

Here, as in (3.3-3.4), we need to identify the UV cut-offs with the constant of integration xq. In order to determine 
the fundamental mass scale parameter ft, which characterizes the region where confinement, DCSB and other non- 
perturbative effects are dominant, we propose to use the following bounds for the pion decay constant in the chiral 
limit 

87.2 < F° < 93.3 (MeV). (4.4) 

The pion decay constant in the chiral limit, evidently, can not exeed its experimental value, so the upper bound in 
(4.4) is uniquely well-fixed. The lower bound in (4.4) is fixed from the chiral value of the pion decay constant as 
obtained in Ref. [42], namely F° = (88.3 ± 1.1) MeV, which, obviously , satisfies (4.4). The value F n = 92.42 MeV, 
advocated in Ref. [43], also satisfies these bounds. We think that chosen interval covers all the realistic values of the 
chiral pion decay constant. In any case it is always possible to change lower bound to cover any requested value of 
the pion decay constant in the chiral limit. This bound is chosen as unique input data in our numerical investigation 
of chiral QCD. The pion decay constant is a good experimental number since it is a directly measurable quantity in 
contrast, for example, to the quark condensate or dynamically generated quark mass. For this reason our choice (4.4) 
as input data opens up the possibility of reliably estimating the deviation of the chiral values from their " experimental" 
(empirical) phenomenologically determined values of various physical quantities which can not be directly measured. 
Thus to assign definite values to the physical quantities in the chiral limit is a rather delicate question (that is why 
we prefer to use and obtain bounds for them rather than the definite values). At the same time it is a very important 
theoretical limit which determines the dynamical structure of low-energy QCD. 

What is necessary now is only to determine the constant of integration x - For this aim, let us remind [19,25,26] 
that DCSB at the fundamental quark level can be implemented by the following commutation relation 

{S-\p),l 5 } + = il.m-V 2 ) ± 0, (4.5) 

so that the 75 invariance of the quark propagator is broken and the measure of this breakdown is the double of the 
dynamically generated quark mass function, 2B(—p 2 ). Let us underline that this relation does not depend on gauge 
and renormalization point choices, by definition, i. e . in any gauge and for any renormalization point this relation 
holds. In particular, denoting as usual the dynamically generated quark mass as the inverse of the quark propagator 
(quark self-energy) at zero point, this quantity, 2B(0), can be defined as a scale of DCSB at the fundamental quark 
level, namely 

Acs Bq = 25(0) - 2m d . (4.6) 

It is worthwhile to emphasize that AcsBq and ra d have direct physical sense within our solutions to the quark SD 
equation since our approach to the deep IR (confinement) region is gauge invariant and free of ghost complications 
(that was discussed above in section 2). Thus we analyse our numerical data at a scale where DCSB occurs at the 
fundamental quark level (remember also the discussion in the main body of the text after Eqs. (2.3-2.4)). What 
is needed in this case is only to simply identify our fundamental mass scale parameter ft with KcsBq, i- e. to put 
ft = AcsBq = 2m d . From Eq. (4.2) then it immediately follows that the constant of integration x is equal to 

x = 1.53 (4.7) 

and it is different from the quenched QCD value x™ m (3.11) at which stationary state occurs there. Without the 
existence of the stationary state (localization of the minimum) one has no criterion how to numerically distinguish 
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them. The advantage of this scale-setting scheme is that, on one hand, it is based on exact relation (4.6), on the 
other hand, numerical values of the mass scale parameter p are extracted from the pion decay constant which, as it 
was emphasized above, is a good experimental quantity. Indeed, substituting (4.7) into Eq. (4.1) and combining it 
with (4.4), one immediately obtains values of the fumdamental mass scale parameter p. This makes it possible to 
numerically calculate all physical quantities considered in our work. These results are presented in Table 1 exept for 
the vacuum energy density. 

It is interesting to note that the information on the dimcnsionlcss constant of integration xo is taken from (4.2) 
and (4.6) on account of the identification p = AcsBq — 2m<j, i. e., in fact, inevitable transformation of the pair of 
independent parameters xq, p into the pair of ^csBq{m<i), P takes place within our scale-setting calculation scheme. 
This is also a manifestation of the phenomenon of the "dimensional transmutation" [44], which occurs whenever a 
massless theory aquires masses dynamically. It is a general feature of spontaneous symmetry breaking in field theories. 

V. NUMERICAL RESULTS AND CONCLUSIONS 

That the QCD ground state (vacuum) has a stationary state for pure gluodynamics (quenched QCD) is a direct 
consequence of the existence and importance of the nonperturbative q~ 4 -type quantum, topologically nontrivial dy- 
namical excitations of the gluon field configurations. Moreover, the inclusion of quark degrees of freedom lifts the 
nonperturbative YM vacuum from its stationarity. Our values for the chiral QCD vacuum energy density are: 

e = -(0.01425 - 0.00196Y/) GeV 4 , (5.1) 
e = -(0.01087 - 0.0015OZV/) GeV 4 , (5.2) 

where the first and second values are due to the upper and the lower bounds in (4.4), respectively. As was mentioned 
in the Introduction, there exists also a contribution to the vacuum energy density at the classical level given by the 
instanton-type nonperturbative fluctuations of gluon fields. Within RILM [5] for dilute ensemble it is: 

ej = - -n = -- x 1.0 fm- A = -(0.00417 - 0.000257V/) GeV 4 , (5.3) 

where 6=11 — (2/3) Y/ is the first coefficient of the /3-function and n is the density of the instanton-type fluctuations 
in the QCD vacuum. This expression was postulated via the trace anomaly relation using the weak coupling limit 
solution to the above mentioned /3-function as well as the phenomenological value of the gluon condensate [45]. 
Recently, in quenched (Nf = 0) lattice QCD by using the so-called "cooling" method the role of instantons in the 
QCD vacuum was investigated [46]. In particular, it was found that the instanton density should be n = (1 + 6)fm~ 4 , 
where 5 ~ (0.3 — 0.6), depending on cooling steps. We find that our values for the vacuum energy density (5.1-5.2) are 
an order of magnitude bigger than the instanton-based RILM (in various modifications) [5,46] can provide at all, Eq. 
(5.3). This can be easly understood indeed since our approach is relevant in the strong coupling limit while instantons 
are weak coupling limit phenomena. 

Stationarity of the nonperturbative YM vacuum and the numerical results for the chiral QCD vacuum energy 
density (5.1-5.2) which for the first time have been calculated from first principles (not postulated via the trace 
anomaly relation as (5.3)) are our main results.]^ It is also clearly shown in our work that precisely the nonperturbative 
q~ 4 -type quantum, topologically nontrivial dynamical excitaions of the gluon field configurations in the QCD vacuum 
is mainly responsible for quark confinement and other nonperturbative phenomena (at least at the fundamental quark 
level). This is our answer to the question arised by DeGrand et al. in Ref. [12] and by Negele in Ref. [46]. At the 
same time, it is very interesting to directly identify this type of field configurations by lattice simulations. 

Recently the effects of nonperturbative QCD in the nucleon structure functions were investigated. A universal 
mass scale parameter m a ~ 470 MeV of the nonperturbative dynamics in QCD was obtained [47]. This is in rather 
good agreement with our numerical results for the physical mass scale parameter p (see Table 1) taking into account 
completely different physical observables have been analysed. The chiral symmetry breaking scale within exact 



7 We use the terms quenched QCD, YM theory and pure gluodynamics as synonyms, as well as the ground state and the 
nonperturbative vacuum. However the ground state of the nonperturbative YM vacuum means its stationary state within our 
terminology. The existence of the stationary state in the nonchiral (real) QCD vacuum due to g _4 -type excitations of gluon 
fields remains an open question yet. 







renormalization group approach [48] is found to be around (400 — 500) MeV. Again rather close to our numerical 
results for p. 

Recent quenched lattice result for the chiral condensate which does not depend neither on the scale nor on the 
renormalization scheme [49] (in our notations), (qq) Q = —(206 MeV) 3 , nicely satisfies our bounds (see Table 1). For 
simplicity, here and below we show only central values, omitting errors due to statistics, renormalization and lattice cal- 
ibration [49] . Recalculated at the conventional scale 2 GeV (with the help of expressions given in [49] for MS scheme) 
it becomes (qq) (2 GeV) — —(245 MeV) 3 . This is in fair agreement with recent phenomenological determination from 
QCD sum rules [50], (qq) (l GeV) = -(229 MeV) 3 , which at the scale 2 GeV becomes (qq) ( 2 GeV ) = ~( 242 MeV) 3 . 
Thus our bounds, formally recalculated at the scale 2 GeV, (230 MeV) 3 < -{qq) (2 GeV) < (245 MeV) 3 , being in 
fair agreement with the above discussed phenomenological and lattice data, favor a large quark condensate advocated 
by the chiral perturbation theory (see discussion in [51]). 

In Table 1, the nonperturbative effective scale for quenched QCD, defined as p§ = x™ ln p 2 , is denoted as A^ while 

for the full QCD, defined as p 2 — xop 2 (and numerical value of xq is given in (4.7)), is denoted as A^JJp . As a subject 
for brief discussion, let us note that an alternative interpretation of these scales is also possible. It is obvious that p 
can be identified with the scale at which chiral symmetry is broken dynamically at the fundamental quark level, i. e. 
p = AcsBq- This was the heart of our setting-scale calculation scheme. It makes sense in addition to identify our 

effective scales with confinement scales, i. e. to put A^ = A[? = Ay M and A^j' = A^ V/ ' = Aq CD for quenched 
(YM) and full QCD, respectively. If QCD confines color, then these two scales in principle need not to be the same, 
while p = AcsBq remains unique as coming from full QCD. Thus within this interpretation there is a hierarchy of 
scales at the fundamental quark level, namely Ay M > Aq CD > AcsBq- The ratio between them are given by the 
corresponding constants of integration. At the same time all these scales are less than the chiral symmetry breaking 
scale at the hadronic level, A x sb < 47rF ff ~ 1.2 GeV [52] (see also Ref. [53] and references therein). At this stage, 
the relation of our scales to the QCD scale parameter, Aqcd, remains, of course, unanswered. We hope to shed some 
light on this possible relation taking into account UV terms in the quark SD equation (2.3-2.4) where Aqcd naturally 
appears via the perturbative renormalization group logarithms. 

For the sake of completeness, let us present the numerical values of basic integrals (containing explicitly quark 
degrees of freedom) which allows one to check our numerical results. These integrals are: J g (1.53, 0) = 0.46, which is 
given in (3.4), 7(1.53, 0) = 0.6876, which is given in (4.1) and J co „(1.53, 0) = 0.6837, which is given in (4.3). Evidently, 
these integrals can be calculated to any requested accuracy. And finally, the numerical values of the bag constant 
B, defined as the difference between the perturbative and nonperturbative vacua [54], are given now by the relation 
B = — e (since the perturbative vacuum is normalized to zero) and can be explicitly evaluated from (5.1), 



and from (5.2), 



B(N f = 0) 


= 0.01425 GeV* = 


- (345.5 MeV) 4 


= 1.85 GeV/fm 3 , 


(5.4) 


B(N f = 


1) = 0.0123 GeV 4 


= (333 MeV) 4 


= 1.60 GeV/fm 3 , 


(5.5) 


B(N f = 2\ 


) = 0.0103 GeV 4 = 


(318.6 MeV) 4 


= 1.34 GeV/fm 3 , 


(5.6) 


B(N f = 3) 


= 0.00838 GeV 4 = 


(302.5 MeV) 4 


= 1.09 GeV/fm 3 , 


(5.7) 


B(N f = 0) = 


= 0.01087 GeV 4 = 


(322.9 MeV) 4 - 


= 1.40 GeV/fm 3 , 


(5.8) 


B(N f = 1) = 


= 0.00937 GeV 4 = 


(311.1 MeV) 4 : 


= 1.22 GeV/fm 3 , 


(5.9) 


B(N f = 2) -- 


= 0.00787 GeV 4 = 


(297.8 MeV) 4 ■■ 


= 1.02 GeV/fm 3 , 


(5.10) 


B(N f = 3) 


= 0.00637 GeV 4 = 


(285.5 MeV) 4 


= 0.83 GeV/fm 3 . 


(5.11) 



The numerical investigation of chiral QCD topology within our model and its comparison with instanton based 
models is the subject for subsequent paper. 
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FIG. 1. The effective potential £l g in (3.10) as a function of xq. A local minimum appears at x™ ln — 2.20 and zero at 
x 1 q = 0.725. Qg falls as ~ —Xq 1 asymptotically at large xq. 
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Table 1: Basic chiral QCD quantities in MeV units 



F° 

IT 


n 


m d 


(qq)o 3 


A (o) 


1V eff 


87.2 


540.6 


270.3 


-201.7 


802 


668.7 


93.3 


578.4 


289.2 


-215.8 


858 


715.5 
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